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Abstract 

We consider the dynamics of N boson systems interacting through a pair potential N~ 1 V a (a?j — 
Xj) where V a (x) — a~ 3 V(x/a). We denote the solution to the ./V-particle Schrodinger equation by 
ipN,t- Recall that the Gross-Pitaevskii (GP) equation is a nonlinear Schrodinger equation and the 
GP hierarchy is an infinite BBGKY hierarchy of equations so that if u t solves the GP equation, 
then the family of fc-particle density matrices {<8>fcWt, k > 1} solves the GP hierarchy. Under the 
assumption that a = N~ e for < e < 3/5, we prove that as N — > oo the limit points of the 
fc-particle density matrices of ipN,t are solutions of the GP hierarchy with the coupling constant 
in the nonlinear term of the GP equation given by J V(x)dx. The uniqueness of the solutions to 
this hierarchy remains an open question. 

1 Introduction 



Consider N bosons in a three dimensional cube A with the periodic boundary condition and volume 
one. The bosons interact via a two body potential 

V a (x) = \v(x/a). 
cr 
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We assume the potential V to be smooth, positive, and with compact support. The parameter a 
determines the range and the strength of the potential: a and iV will be coupled so that a — > as 
N — > oo. Thus the potential V a converges to a Dirac 5-function. The A-body Hamiltonian for the 
N weakly coupled bosons is thus given by 

N 1 

H N = - Aj + - V a ( Xl - x 3 ). (1.1) 

The density of the bosons in the cube, p, is clearly equal to N. Any given particle typically 
interacts via the potential V a with a 3 N other particles. If a S> A~ 1//3 , there are a lot of interactions 
among the N bosons. In this case, the ground state of the system contains, to the leading order in 
N, no correlation among the particles. Thus the ground state wave function is a product function to 
the leading order. In fact, the same conclusion holds as long as a > A -1 . The correlations among 
the particles become important in the leading order only when a ~ A -1 . This choice corresponds to 
the so called Gross-Pitaevskii scaling limit, as pointed out by Lieb, Seiringer and Yngvason [B] (see 
[7j for a review). For this choice of scaling, i.e., a ~ -/V -1 , we study the dynamics of the Bose gas in 
[I]. In this paper, we consider the cases a = N~ £ for < e < 3/5. 

In the following we denote by x a general variable in the box A. On the other hand x = 
(x\, . . . xn) denotes a point in A^. We will also use the notations x^ = (xi,...Xk) £ A fc and 
X7V-/c = (scjfe+i, • • • x N ) G A N ~ k . 

The dynamics of the Bose system is governed by the A^-body Schrodinger Equation 

id t ipN,t = H N ti) N>t ■ (1.2) 

Here the wave function ipN,t £ -^(A^), the subspace of L 2 (A N ) consisting functions symmetric with 
respect to permutations of the N particles. More generally we can describe the N body system by 
its density matrix jN,t- A density matrix is a positive self-adjoint operator 7 acting on L^(A N ), with 
TV7 = 1. The density matrix corresponding to the wave function ipN,t is given by the orthogonal 
projection onto tpN,t, i-e., ^N,t = \ipN,t){ipN,t\- Quantum mechanical states described by orthogonal 
projections are called pure states. In general a density matrix is a weighted average of orthogonal 
projection (mixed states). The time evolution of the density matrix 7^ is given by 

id t lN,t = [HN,7N,t] (1-3) 

which is equivalent to the Schrodinger Equation (|1.2|) . 

It is useful to introduce the marginal distributions corresponding to the density matrix 7jvt- For 

(k) 

k = 1, . . . , N, the /c-particle marginal ^ N t is defined through its kernel by 

7^1( x fc; x fc)= / dx Af _ A .7 A r it (x fc ,x Ar _ fc ;x / fc ;x iV _ fc ) (1.4) 
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where xl = (x^, . . . ,x' k ) and where 7jv,t(x;x') denotes the kernel of the density matrix jN,t- From 
TrTiV,* = 1, it immediately follows that 

Try™ = 1 (1.5) 

for every k = 1, . . . ,N. 

Using (|1.3j) (and the symmetry of ^N,t with respect to permutations of the N particles) we find 
that the evolution of the marginal distributions of jN,t is determined by the following hierarchy of 
iV equations, commonly called the BBGKY Hierarchy: 

k 



3=1 

1 * 

+ n - x *) - V M - 4))7S,l(x fc ;x' fe ) (1.6) 

+ (- 1 ~ E J dx k+i{V a {xj - x fc+ i) - V a {x'j - x fc+1 ))7^ 1) (x fc ,x fc+ i;x / fc ,x fc+1 ). 



3=1 

Rewriting this hierarchy in integral form yields 

k 



1n\ ( x ^ ; 4 ) = 7/?,o ( x fc ; 4 ) - iV / ds (- a Xj + a x , ) 7 ^ ( Xfc ; 4 ) 

i=l Jo 

• to rt 

N m h 

J ds J dx k+i {Va{xj - x k+l ) - V a (x' j - x A . +1 ))7^ 1) (x fe ,x fc+ i;x / A .,x fc+ i). 

(1.7) 

Letting a = iV~ e and considering the limit N — > 00, the BBGKY Hierarchy converges formally to 
the following infinite hierarchy of equations: 

k . 

7 f ) (x fc ;x' fc )= 7 f (x fc ;x' fc )-z^ f ds (-A Xj + A x , ) 7 i fc) (x fc ; x' fc ) 

i=i 

-i&Y] / ds / da; fc+ i(5(xj - x fc+ i) - 5{x' j - x k+1 ))^ s k+1 \x k ,x k+1 ;x' k ,x k+1 ), 



(1.8) 
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where we defined 



b = / dxV(x) 



(recall that V a (x) = (l/a 3 )V(x/a)). We will call (|1.8|) the infinite BBGKY hierarchy, or the Gross- 
Pitaevskii (GP) hierarchy. Note that (|1.8j) has a factorized solution. The family of marginal distri- 
butions 7 t (xjt; x' fc ) = nj=i 4>t{xj)4>t{x'j) is a solution of ()1.8|) if and only if the function 4> t satisfies 
the non-linear Schrodinger Equation 

idtMv) = -&Mv) + b\Mv)\ 2 My) • (1-9) 

This is the Gross-Pitaevskii (GP) equation, except that the coupling constant in front of the non- 
linear interaction is given by b. In the standard GP equation, the coupling constant is 87rao where 

00 is the scattering length of the unsealed potential V{x). 

The aim of this paper is to prove the convergence of solutions of Q1.6|) to solutions of (|1.8[) : 
more precisely we will prove that the sequence Tnj = {7yvi}fcLi nas & t least one limit point = 

{7oo\}fc>i w ith respect to some weak topology, and that any weak limit point satisfies the 

infinite hierarchy (|1.8|) . For dimension d = 1, the convergence to the GP hierarchy (|1.8|) for the delta 
potential was established by Adami, Bardos, Golse and Teta in 1 . If the potential has a weaker 
singularity, the convergence to the infinite BBGKY hierarchy was proved in [2]. 

In [I] we proved the convergence to the GP hierarchy with the coupling constant 87rao when 
a ~ A -1 . Together with the result of the present paper, this means that the coupling constant 
changes from 87rao to b when a changes from A -1 to N~ £ , for e < 3/5. In fact, the coupling constant 
should always be given by b, as long as a = N~ £ , with e < 1. This fact follows from the observation 
that the potential N~ 1 a~ 3 V(x/a) has scattering length of order 1/A for every choice of a = N~ e , 
e > 0. Hence the ground state correlations live on the scale 1/iV, for all e > 0. As long as e < 1, 
these correlations cannot affect the coupling constant because the potential live in the much larger 
scale a = N~ £ 3> A^ 1 . In other words, since the two-particles correlation function is of the form 

1 — const/(A|x|), for |x| > A 7 ' -1 , we have 

which equals b in the limit N — » oo, if a ^> N . Thus, a ~ A -1 is the only scaling for which the 
coupling constant is given by the scattering length of the unsealed potential V. 

We remark that in [3], we need to use a modified version of the Hamiltonian (|1.1|) . in which 
two-body interactions are removed, whenever too many particles come into a small region. In the 
present paper, this assumption is not needed. Moreover, we prove a much stronger a-priori estimate 
on the limiting density matrices, see (|2.16|) . 

Since the 5-function cannot be bounded by the Laplace operator, we are unable to prove the 
uniqueness of the solutions of ()1.8|) . Hence we cannot conclude the propagation of chaos. The best 
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known result in this direction is [3] , which covers the case of a Coulomb singularity in the potential in 
d = 3. Previously, uniqueness was proved by Hepp |Hj and Spohn [3] for bounded potential; Ginibre 
and Velo _5 had a completely different approach for quasifree states. 



2 The Main Result 

(k) 

Since our main result states properties of limit points of the sequence ^y N t for N — > oo, in order to 
formulate it, we need to specify a topology on the space of density matrices. 

Quantum mechanical states of a &-boson system can be described by a density matrices 7^: 7W 
is a positive, trace class operator, with trace normalized to one. We can also identify 7W with its 
kernel and consider it as a distribution in L 2 (A k x A fc ): in fact, since 7^ is a positive operator with 
trace equal to one, its Hilbert Schmidt norm is also bounded by one, and 

|| 7 {fc) ||2 := J dx fc dx' fc | 7 (fc) (x fc ;xi.)| 2 < 1. (2.11) 

For T = {7 (fc) } fc >i G ©fc>iL 2 (A fc x A fc ) we define the norm 

00 

||r|| H _ :=^2- fc || 7 ( fe )|| 2 , (2.12) 
k=i 

and we put 

H_ := {r € 0L 2 (A fc x A fc ) : ||r|| H _ < 00}. 
k>l 

Analogously, we define 



H + = {T (£ ff)L 2 (A k x A fc ) : lim 2 fc ||7( fe )|| 2 = 0}, 

k>l 

and we equip H + with the norm 

||r|| H+ =su P 2 fc ||7( fc )|| 2 . 

k>l 

The Banach space J3_ is the dual space to H + . A sequence Fn = {7/^^}fc>i G H- converges to 
Poo = {7oo^}fc>i G H- in the weak* topology if and only if 

J^E / dx fc dxU (fc) (x fc ,x' fc )( 7 {f ) (x fc ,x' fc )-7W(x fc ,x' A; )) =0 (2.13) 

for all J = {J( fe )}fc>i G H + (this is actually equivalent to convergence for each fixed k). We will 
denote by C([0, T], H-) the space of functions of t G [0, T] with values in i2_ which are continuous 
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with respect to the weak star topology on H-. Since the space H + is separable, we can fix a dense 
countable subset in the unit ball of H + , denoted by {Jj}j>l- Define the metric on H_ by 

oo oo „ 

p(T, f ) := E 2 "1 E / dx * dx U (fc) 4) b (k) - ¥ k) (x fc ; x> fc )] | . (2.14) 



i=l fc=i 

Then the topology induced by p(., .) and the weak* topology are equivalent on the unit ball B- of 
H_. We equip C([0,T],H_) with the metric 

p(r,f):= sup p(T(t),f(t)) . (2.15) 

0<t<T 

We are now ready to state our main theorem. 

Theorem 2.1. Assume the potential V(x) is positive, smooth, and has compact support, and set 
V a (x) = a~ 3 V(x/a). Suppose that a = N~ £ , for some < e < 3/5. Choose an initial density matrix 
7n,o such that 

TrH k NlNfl < C k N k 

for some constant C and for all k > 1. Let rjy,o = {7jvo}fc=i ^ e ^ e / a ™/y °f marginal distributions 

corresponding to the initial density matrix 7jv,o- Fix now T > and denote by T^,t = {7w1}j f or 
t G [0, T], the solution to the BBGKY Hierarchy \l.b}) corresponding to the initial data Ynq. 

i) The sequence Yjsit *s compact in C([0, T],i2_) with respect to the metric p. 

ii) Let Tqo^ = {7^\}fc>i G C([0, T], 6e any Zwnii pom£ o/Tjv,t respect to the metric p. 
Then there is a constant C such that 

Tr\S 1 ...S kl t]tSk-..S 1 \<C k (2.16) 

/or every k > 1. -f/ere we use i/ie notation Sj = (1 — A.,) 1 ' 2 . 
Too^ is non trivial. In particular we have 

Trit]t = 1 (2-17) 

/or every t G [0, T] and for every k > 1. 

raj Assume h r (x) = r~^h{x/r) for any h G Cq°(A) wii/i f^h = 1. Then, for any k > 1 and 
£ G [0, T] ; i/ie Zimit 

lim / dxj b+1 dx fe +i/i r (a;' A . +1 - ajfc+iJJvCafc+i - Xj)^\(x k ,x k+1 ;x' k ,x' k+1 ) 
ry^oj (2.18) 

_ ( fc ) 1 1 \ 

— ■ 7oo,t l x fc ) x i > x fc 1 ^ j J 
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exists in the weak W 1 ' 1 (A fe x A k )-sense and defines 7™(xi, Xj\ x' k , Xj) as a distribution of 2k 
variables . 

v) ^oo,t satisfies the infinite Gross-Pitaevski Hierarchy U.8\) in the following sense: For any 
jW(x fc ,x' fc ) e W 2 >°°(A fc x A k ) we have 

J dx k dx' fc J (fc) (x fe ,x^.)7^ t (x fc ,x / fc ) = j dx k dx' k J {k \x k ,x k )^] Q (x k ,x' k ) 




(2.19) 

Here the action of the 5 -functions on 'y( fc+1 ) is well defined (through a regularization of the 
5 -function) by part iv). 

3 Energy Estimates 

The main tool in the proof of Theorem 12.11 is the following proposition, which proves bounds for 
the L 2 -norm of the derivatives of a wave function ip in terms of the expectation of powers of the 
Hamiltonian Hn (defined in (jl.ljO in the state described by ip. The proof of this proposition requires 
some standard Sobolev-type inequalities, which are collected, for completeness, in Appendix 1X1 

Proposition 3.1. Suppose V is smooth and positive. Put V a (x) = a~ 3 V(x/a) and assume a = N~ s , 
with < e < 3/5. Put 

N 1 

j=l e^m 

Fix fceM and < C < 1. Then there is N = N (k, C) such that 

(V, HU) > C k N k (V, SlS 2 2 . . . Sfa) (3.20) 

for all N > Nq and all ip € D(H^) (ip is assumed to be symmetric with respect to any permutation 
of all its variables). 

1 Here and in the following W p ' q (A k x A fc ) denotes the usual Sobolev space over A k x A k . 
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Proof. The proof of the proposition is by a two step induction over k. For k = and k = 1 the claim 
is trivial (because of the positivity of the potential). Now we assume the proposition is true for all 
k < n, and we prove it for k = n + 2. To this end we apply the induction assumption and we find, 
for N > N (n,C), 

(tl>, H n + 2 ^) = ty, H N H n N H N ^) > C n N n ^H N Sf . . . S 2 n H N ^) . (3.21) 

We put 

n N 
3=1 j<m 

with Vj m = a~ 3 V((xj — x m )/a). Then we have 

{i,H N Sl...S 2 n H N i,)= J2 (^,S 2 1 S 2 ...S 2 n S 2 ^) + J2 (^ 1 S 2 S 2 ...S 2 n H^^) + c.c.) 

ji,32>n+l j>n+l 

+ (4j,h^sI...sIh^4j). 

Since H^S 2 . . . S 2 H^ > 0, we find, using the symmetry with respect to permutations, 

m N S 2 . . . SlH N ^) >(N-n){N-n- S 2 . . . S 2 n+2 ^) + (2n + 1)(N - n)ty, SfS 2 . . . S 2 n+1 ^) 

n(n + 1)(N — n) , ^ 2 ,, N 

+ 2iV ^' 12 1 ' ' ' n+l ^ + CX ^ 

(n + l)(iV-n)(iV-n- 1) , Tr _ 2 n2 M N 

+ " ^ ^l,n+25f • . . S 2 n+l i>) + C.C.) , 

(3.22) 

where c.c. denotes the complex conjugate. Here we also used that (ip, Vj m S 2 . . . S 2 +1 ip) > if 
j,m > n + 1, because of the positivity of the potential. Next we consider the term on the second 
line of (|3.22|) : note that this term vanishes if n = 0, so we can assume n > 1. Then we have 

(V>, V 12 Sf . . . Sl +l tp) + c.c. = (V>, S n +i • • • S 3 Vi2S 2 S 2 2 S 3 . . . S n+1 ip) + c.c. 
= (V, S n+1 . . . S 3 Vu(l + p?)(l + p 2 2 )S 3 . . . S„+iV) + cc. 

> 2(tp, S n+ i . . . S 3 Vi2plS 3 . . . Sn+iip) + (ip, S n+ i . . . S 3 Vi2pIp%S 3 . . . S n+ itp) + c.c. 

> 2(ip, s n+ i . . . s 3 w 12 p 2 s 3 . . . Sn+iVO + (i 3 , s n +i ■ ■ ■ s 3 p 2 Wi 2 pip 2 s 3 . . . s n+ i4>) 

+ (tj>, S n+ i . . . S 3 W 12 p 2 p 2 S 3 . . . S n+1 tp) + c.c. 
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where VV12 = a 4 (W)((xi — X2)/a). Applying Schwarz inequality we get 
{^VxiSl . . . S 2 n+l ^) + c.c. 

> - 2{aity, S n+1 ...S 3 \VV 12 \S 3 ... 5 n+ iV) + a^ty, S n+1 . . . S 3 \p 2 \\VV 12 \\p 2 \S 3 . . . S„+iV)} 

- {a 2 (ip, S n+1 . . . .S3IP2II W12IIP2IS3 • • • S n+1 ip) 

+ a2 1 ('i/f,iS r n+ i . . . S3 |j9 2 1 Ipi 1 1 VV12 1 |pi I |P2 I ^3 • • • S n+ iip)} 

- {a 3 (ip,S n+1 . . . S3IW12IS3 • • • S n+1 ip) + a^ 1 (ip,S n+1 . . . S 3 \p 2 \pj\W 12 \pl\p 2 \S 3 . . . S n+1 ip)} . 
Using Lemma lA. II we find 

&,V 12 S 2 . . . S^) + c.c. 

> - C{a l a" 1 (i), S n+ i . . . S 3 SfS 2 S 3 . . . S n+1 ip) + cq l a~ 2 (ip, S n+ i . . . S 3 Sfp^S 3 . . . S n+1 ip)} 

- C{a 2 a~ 2 (ij), S n+ i . . . S 3 SfplS 3 . . . S n+ iip) + a 2 1 a~ 2 (ip, S n+ i . . . S 3 S\p\plS 3 . . . S n+l ip)} 

- C{a 3 a~ 1 ('ip, S n+ i . . . S 3 SlSlS 3 . . . S n+1 ip) + aJ 1 a" 4 (V>, S n+ i . . . S 3 p\p\S 3 . . . S n+ iip)} 

> _ CN- 3 / 2 a- 5 / 2 {N 2 (iP, S 2 S 2 2 S 2 3 . . . S 2 n+1 iP) + Nty, Sfs! . . . S 2 n+1 ^)} 

(3.23) 

where we optimized the choice of aci, a 2 , a 3 . As for the last term on the r.h.s. of 1)3.22(1 we have 

O, Vi, n+2 S 2 . . . S 2 +1 ifj) + c.c. = (ip, S n+1 . . . S 2 V^ n+2 SlS 2 . . . S n+l tp) + c.c. 

> (if), S n+ i . . . S 2 Vi >n+2 p\S 2 . . . S n+ iifj) + c.c. 

> (ip, S n+ i . . . S 2 VVi, n+2 piS 2 . . . S n+ iip) + c.c. 

> - a(ip, S n+1 . . . S 2 \Wi, n+2 \S 2 . . . S n+ iip) 

- a" 1 ^,^! . . . 5 2 |pi||VVi in+ 2||pi|S'2 . . . S n+1 ip) 

> - C{aa~ l + a~V 2 ) (i/>, S 2 . . . S 2 n+2 ^) 

> -Ca- 3 / 2 (iP,S 2 ...S 2 n+2 iP). 

Inserting last equation and (|3.23|) in the r.h.s. of (|3.22|) we get 

m N S 2 . . . S 2 H n tP) > (N - n)(N - n - 1) (l - ^ - ^5^) Sf... S 2 n+2 ^) 

+ (2n + 1)(N- n) (l - j^j^j (A SfS 2 . . . S 2 n+l ^) . 

Since a = N~ £ with e < 3/5, we have A r3 / 2 a 5//2 S> 1 and Na?/ 2 ^> 1. From the last equation, for 
any fixed C < 1 and n £ N, we can find A^o so that 

(V, H N S 2 . . . S 2 n H N ^) > C 2 N 2 (t(;, S 2 . . . S 2 n+2 ^) (3.24) 

for every A^ > Nq. This, together with (|3.21j) completes the proof of the proposition. □ 
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Corollary 3.2. Suppose the initial density matrix 7jv,o satisfies 

TrH k NlNfi < C k N k . 

Let 7jv,t be the solution of (|1.3|) . and let {r/$t^k=a ^ e ^ e corresponding marginal distributions. Th 
for any C > C\ and any k € N there is Nq = No(k, C) such that 

for allteR and all N > N . 



en, 



TrS 1 . . . S k ^ t S k ...S 1 <C k (3.25) 



4 Proof of the Main Theorem 

In this section we prove our main result, Theorem 12, 11 To this end we will make use of the following 
lemma. We use here the notation 

<JW 7 (fc) > = / dx fe dx' fc j( fc )(x fc ,x' fe )7 (fe) (x fc ,xl,) . (4.26) 

Lemma 4.1. Fix k > 1 and G W 1>00 (A k x A k ). For (3 > set 5p{x) = (4tt/3 3 /3)^ 1 x(|^| < &)■ 
Then we have 

k . 

(J (k) ^t) = (^.^-''EjJ ds / dx,dx' fe J«(x fe ,xj t )(-A :c . + A 2 ,)7^( x *>4) 



3=1 



J ds J dx fc dx fc dx fe+i j(fc) ( x ^' x fc)^/3( x i ~ ~ ~ ( 4 27) 



+ t (kO(j3V 2 ) + kO(a 1 / 2 ) + k 2 (j^j^j ) sup Tr|5i 527^52^1 . 
Proof We start from the BBGKY Hierarchy (fT77|) . After multiplying with j( fc ) E L 2 (A fc x A fc ) we 
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get 



k ^ 

(J^^nI) = (^,7g)-'E / ds I dx fc dxU (fc) (x fc ,x' fc )(-A x , + A x ,.hS(x fc ,x' fc ) 

j=1 J° J 



l 

N 



o 



^2 I ds I dx fc dx' fc j( fc )(x fc ,x' fc )(y a (x i -x e ) -Vaix'j - x' e ))^ s (x k ,x' h 



k 



" 1 ~ J ds J dxkdx' k dx k +i J (fc) (x fc ,x' fc ) (V a (iCj - a^+l) - - x k+1 )) 



3=1 



x T Ns '{x k ,x k+ i,x k ,x k+ i). 



(4.28) 



Next we estimate 

J dx fe dx^ jWfa^Vaixj -x/)75fi( x ki4)| = ^^^(xj -x £ ) 7 $i 

< ll-sj 1 jC*)^- 1| UsrVafo - x^sj 1 1| 1| TrSjSr/filSjSt ( 4 - 29 ) 

<C k a- 1 / 2 TrS 1 S 2 'y%lS2Si 

where we used that, bv LemmaO ||(V a ) 1/2 (^-^)5' j " 1 || < Ca^ 1 / 2 , and 1 1 (K ) 1/2 j - ^ ) 1 S7 1 1 1 < 
C. Moreover we used that WS^J^SjW < \\J^Sj\\ and that 

||j( fe )^-|| 2 < ||j( fc )S 2 (j( fc ))*|| < Tr J^S](jWy = J dx fc dx' fc (|j( fc )(x fc ,x' fc )| 2 + |V,jW(x fc ,x' fc )| 2 ) 

(4.30) 

which is bounded, because of the finiteness of the volume, and because, by assumption G 
W 1 '°°(A fc x A k ). 

In the same way we can bound the contribution arising from the term V a (x'j — x'g), and so we find 



J ds / dxfedx fc J {k) {^k^ k ){V a {xj -x e ) -Vaix'j -x^))7^ ) s (x A; ,x / fc ) 



< tttJo SU P TrSi^Twl^Si ■ ( 4 - 31 ) 
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Analogously we also get 

k „ t 



k 
N 



y2 ds dx fc dx' fc dx fc+ i J W (x fe , x' fe ) (V a (xj - x k+x ) - V a (x'j - x k+1 )) 
3=l Jo J 



(fc+l)/ / X 



< -^WJ^WWSkhVaixj - Xfc+i)5 fe ^5- ^IIIS- *|| sup Tr 5i5 2 7^gS , 2'S , i 
7V ' se[o,t] 

Cfci 



(4.32) 



< 



sup Tr5i5 2 7^S 2 Si 
Na 1/2 aG [ ,t] 



Applying Lemma 14 . 21 twice (once with /5 2 = a and once with /3 2 = /?; in both cases with /3j =0), we 
have 



" - s [x-k,x k+ i,x k ,x k+1> 



J dxfedx^. J (fc) (x fe ,x' fc ) (V a (xj -x k+ i) -b5p(xj -x k+1 ))j^ : 

<C{a 1 ' 2 + ^! 2 )^S 1 S 2l ^] s S 2 S 1 

for some constant C which only depends on j( k \ but is independent of N, of 0, and of s S [0, t]. □ 

The following lemma is used to regularize the action of the <5-function. It was already used in the 
proof of Lemma l4,ll It will be used again to prove the convergence to the infinite BBGKY hierarchy. 
Its proof can be found in 0j (see Proposition 8.1). 

Lemma 4.2. Suppose 8a(x) is a radially symmetric function, with < 8g(x) < C(3~ 5 x(\ x \ < P) 
and J 5p(x)dx = 1 (for example 8a(x) = 0~ 3 h(x/ (3), for a radially symmetric probability density 
h{x) supported in {x : \x\ < 1}). Then, for any 

j(k) e w i,°o( A k x A fc) and for any smooth function 
^/( fe + 1 )( Xfc+1 ; x' fc+1 ) corresponding to a (k + V) -particle density matrix, we have, for any fixed j < k, 

J dx k dx' k dx k+1 dx k+1 J (fc) (x fc ;x' fc ) {5 Pl {x' k+1 - Xk+i)b~p 2 {xj - x k+x ) - 5(x k+1 - x k+ i)S(xj - x k+1 )) 



X7 



(fc+i) 



( X fe; X k+ l,TL k , X k+l ) 



< CipWco + HVjJlloo] (Jh + Vfa) TrlSjSk+a^SjSk+il . 



(4.33) 



4.1 Compactness of the sequence rV(i) 

The aim of this section is to prove part i) of Theorem 12.11 
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Proof of Theorem \2.1\ part i). First of all we note that the sequence IV,* = {7^+}feLi i s uniformly 

(k) (k) 

bounded in H-. In fact, from j Nt > and Trj^ri = 1 it follows immediately that 



ir 



N,t\\H- 



^2- fc f dx k dx' k \^ t (x k - X ' k )\ 2 <l. 
k>l 



(4.34) 



Next we prove that the sequence r^t is equicontinuous in time with respect to the metric p (see 
(|2.14|) ) defined on H_. To check equicontinuity we use the following lemma, whose proof can be 
found in [3] (see Lemma 9.2). 

Lemma 4.3. The sequence Tnj = {7yvt}fcLi> N = 1,2, . . . satisfying \4-^4\) ^ s equicontinuous on H_ 
with respect to the metric p if and only if for every fixed k > 1, for arbitrary £ W 1 '°°(A k x A fc ) 
and for every e > there exists a 5 > such that 



< e 



(4.35) 



whenever \t — s\ < 5. 

So, in order to prove that IV t is equicontinuous, we choose k > 1, J (fc) G W 1 » 0O (A* x A fc ) and 
e > 0. Then by Lemma 14. 11 we have 

' i=i Js 

+ &E / dr / dx fc d x / fc dx A . +1 J (fc) (x fc ,x' fc )((5 /3 (x i - x fc+ i) - ^(x^ - x k+1 )) 

3=1 Js J 

X 1N,t (Xfc,X fc+ l,X fc ,Xfc+l) 

+ \t-s\ (kOiP 1 / 2 ) + W^ 2 ) + k 2 (j^)) > 



(4.36) 



sup Tr\S 1 S 2 T N ' T S 2 Si\ 

6[0,t] 



Next we note that 



(4.37) 



is uniformly bounded in A r and in r € [s,t], because of Corollary 13.21 and of (|4.3U|) . 

As for the term on the second line of (|4.36j) we note that, for fixed /3 > 0, it is bounded by 

/ dr / dx fc dxfcda; fc+ iJ (A:) (x fe ,x / fc )(^(x i - x k +i) - 8p(x'j - x k+1 ))^^ 1 \y: k ,Xk + i,x k ,x k+1 ) 
<C/3" 3 E / dr||j( fe )||IV 7 ^ 1) <C /9iM(fe) |t- S | 

1 J s 



(4.38) 
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where we used that is a bounded operator (this follows easily from the condition that its kernel 
lies in W 1,00 {K k x A fc ), and from the finiteness of the volume), that the norm of 8a is of order /3 -3 , 

and that the trace of ") Nt is one, for every r and N . From the last three equations we find 



7& 



<C\t-s\ 



(4.39) 



for a constant C, depending on (5,k,J^ k \ but independent of N,t,s. This implies, by Lemma 14.31 
equicontinuity of rjv,t- The equicontinuity of r^t then implies that rV,t is compact in C([0, T],H_) 
by the Arzela-Ascoli Theorem. □ 



4.2 A -priori bounds on ^ 

The aim of this section is to prove part ii) of Theorem 12. II To this end we define a new topology in 
the space of density matrices. 

Denote by £ l (TL) and by K,(H) the space of trace class operators and, respectively, the space 
of compact operators on a Hilbert space 7i. Moreover let H k = L 2 (A k ). For a density matrix 
ry{k) g £ 1 (/f fc ) ) we define the norm 

\\ 1 W\\ Wk =T I \S 1 ...S kl Ws 1 ...S k \ (4.40) 
where Sj = (1 - A^) 1 / 2 . We put 

W fc = { 7 (fe) G ^(fffc) : ||7 (fc) llw fc < oo} . (4.41) 
We consider moreover the space 

A [k) = | T (fc) =Sl ...S k S!...S k : e JC{H k )} (4.42) 
equipped with the norm 

r (fc) iUw = \\s^...s^t^s^...s^\\ (4.43) 

where |.| denotes the operator norm. We have 

M (fc) ,IMU«)* = (w fc ,||.||w fc ). 

The identification of W k as the dual space to implies the existence of a weak star topology on 

Proof of part ii) of Theorem \2.1\ Let = {7^\}fe>i be any limit point of rj\r,t = {l^\} k= i in 
the space C([0,T],H-) with respect to the metric p. By passing to a subsequence we can assume 
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that Tnj ~~ > Tqo^, for N — > oo, w.r.t. the metric p. This implies that, for every fixed t G [0, T] and 
for every k > 1, we have 

7$ - 7S (4-44) 

with respect to the weak topology of L 2 (A k x A fc ). This follows because ||rAr jt ||#_ < 1 and because 
in the unit ball, the metric p is equivalent to the weak* topology of i?_ . Convergence with respect 
to the weak* topology of H- implies then weak convergence in every fc-particle sector L 2 (A k x A k ). 
By Corollary 13.21 there exists a constant C such that 



hjvJllWfc = Tr Si... S k ^ t S k ...St 



(fc) 



< C k 



(4.45) 



for every t G [0,T] and k > 1. By the Banach-Alaoglu Theorem, the sequence 7^ is compact in 
Wfc with respect to the weak* topology. In particular there exists a subsequence Nj — > oo, and 

6 W fc such that 7^ -» 7^ and 



llnSillw* = Tr Si... ShT^tSk ■■■Si 



< C k 



(4.46) 



So, the sequence 7^ t satisfies, for j — > 00, 



(fc) _ 

7jVj,t ~~ > 7oo,t 



7iv. 



^ w.r.t. the weak topology of L 2 (A k x A k ) and 

( fc ) . 



Too t w.r.t. the weak* topology of W k ■ 



(4.47) 



If G L 2 (A k x A fe ) then the operator with kernel given by (which will be still denoted by 
JW) is Hilbert-Schmidt and thus compact: in particular G Ak- Thus, using Q4.47JI . it is easy to 
verify that 

J dx fc dxl. J (fc) (x fc ,x' fc ) 7 2; t (x fc ,x' fc ) = I dx fe dx' fe J( fe )(x fe ,4)7^(x fe ,4) (4.48) 

for every G L 2 (A k x A k ). This implies that 7^ = 7^ as elements of L 2 (A k x A fc ). Thus, from 
l|PS|t . we have 



Tr 



£1 • • • Sk^tSt . ..Si 



< C k (4.49) 

for every i G [0, T] and A; > 1. More precisely one should say that there is a version of 7^ G 
L 2 (A k x A k ) which satisfies this bound (the version we are using here is exactly the density matrix 



□ 
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4.3 Non-triviality of the limit points 

Proof of part Hi) of Theorem \2.1\ Suppose = {j^\}k>i is a limit point of T^t in the space 
C([0, T], H_) with respect to the metric p. From Section f4. 21 we know that, for every fixed t G [0, T] 

(k) / (k) 

and k > 1, 7^ t € V\4 (more precisely there is a version of 7^ lying in the space W&), and that 

(k) (k) 

there is a subsequence Nj — > 00 with 7^ , -» 7^ w.r.t. the weak* topology of This means 
that 

TrJ (fc) (7^-7S)-0 (4.50) 

for j -> 00, and for every j( fc ) G A (recall that j( fc ) G A- if and only if 1 . . . S fc 1 J (fc) S fc 1 . . . 1 
is compact as operator on L 2 (A fc )). Next we note that, because of the finiteness of the volume of A, 
the identity operator is an element of A k (since 5-f 2 . . . S^ 2 is a compact operator on L 2 (A k )), and 
thus 

Tr 7 W = limTr 7 « =1 (4.51) 
because Tr^ t = 1 for every N,t,k. □ 

4.4 Convergence to the infinite BBGKY Hierarchy 

In this section we prove the last two parts of Theorem 12.11 

Proof of part iv ) of Theorem \2.1\ Eq. (|2.18|) follows from part ii) of Theorem 12.11 and from Lemma 
14.21 by the following simple argument. From Lemma 14.21 we find 

J dxfcdxfcdar fc+ i<ia4 + i J (k) (x k ; x' k ) 7^ 1} (x fc , x k+ i; x' fc , x' k+1 ) 

X (6 ri (x' k+1 - Xk+ijSr^Xj - X k +l) - S r2 (x' k+1 - X k+ i)S r ' 2 (xj - Xk+l)) 

< C[\\ JIU + [|V; J||oo ] (ri + r 2 + yjr~[ + Tr \S 3 S k+ll tf ] SjS k+1 \ . 
This implies, by (|2.1H|) . that the sequence 

J dx fe dx' fc dx fe+ idx / fc+1 J (fc) (x fc ;x' fc ) S r (x' k+1 - x k+1 )5 r >(xj - x k+1 ) 

(k+l) ( / 1 x 

X 7oo,i V x fc > x k+l ; X fc , X k+1 j 

has the Cauchy property for r, r' — > and thus converges, if J( fe ) G W 1 '°°{A k x A fc ). □ 



(4.52) 



(4.53) 
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Proof of part v) of Theorem \2.1\ From Lemma 14.11 we find, for an arbitrary G W l,co {A k x A k ) 
and for TV large enough, 

k 



-iby^ ds dx k dx' k dx k+ ij( k \x k ,x' k )(5f3(xj - x k+x ) - 8p(x'j - x k+1 )) (4.54) 

, =1 y 



x 7at,s {^-k,x k+1 ,x k ,x k+l ) 
+ tO{l3 l ' 2 )+to(l) 

where o(l) — > for A" — > oo. By passing to a subsequence we can assume that T^t — * = 
{7oot}fc>l e w.r.t. the metric p. Since ||rjv,t||if_ < 1, and since the metric p on the 

(k) (k) 

unit ball of iT_ is equivalent to the weak* topology, it follows that ^ N t —> 7( ^/ t w.r.t. the weak 

topology of L 2 (A k x A'), for every fixed k > 1 and t G [0,T]. For jW G IT 2 ' 00 (A* x A fc ), we also 
have j( k ' G L 2 (A k x A fc ) (because of the finiteness of the volume). Hence 

<^7^-^S>-° and <^7^o-72il,>-0 (4.55) 

for N — > oo. 

As for the second term on the r.h.s. of (|4.54|) we note that, from G W 2 ' oc (A k x A fc ), it also 
follows that A Xj j( k '(x. k ; x' fc ) and A x / J^(x k ; x' fc ) are elements of L 2 (A k x A fc ). This implies that 

2 / dx fc dx' fc A,J«(x fc ;x' fc ) ^(x^x',) - 7 g s (x fe ;x' fc )) -> (4.56) 
for A" — ► oo, and for every s G [0, t]. By Lebesgue Theorem on the dominated convergence, we find 
fds I dx fe dxU (fe) (x fc ;x;)(-A a ,+A,,)( 7 J ) s (x fc ;x;)- 7 g s (x fc ;x' fc )) ^0 (4.57) 

j=l J J 

for N — > oo and for every fixed t G [0, T]. 

Finally we consider the limit A" — ► oo of the last term on the r.h.s. of l|4.54|) . From Lemma 14.21 
we have 

J dx fc dx' fc dx fc+ i J (fc) (x fe ;x' fc ) (Sp(xj - x k+1 ) - Spty - x k+x )) 7 ^ 1} (x fc , x k+1 ; x' fe , x k+l ) 
= j dx A ,dx' fc dx fc+ idx' fc+1 J (fe) (x fc ;x' fc ) {8p(xj - x k+1 ) - 8p{x' j - x k+1 )) 

x 8 v (x k+ i - 4+i)7vt 1) ( x fc' x fc+i; x fc>4+i) + ( 4 - 58 ) 
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(4.59) 





where 0(rj) is independent of (3,N and s. At this point we can take the limit N — > oo with fixed j3 
and rj. Since G M^ 2 '°°(A fc x A fc ), it is easy to check that, for fixed (3, rj > 0, j( fc )(x fc ; x' fc )(5 / g(xj — 
2^+1)^7(2^+1 — ^fc+i) is an element of L 2 (A fc+1 x A fc+1 ). Hence 

J dx k dx' k dx k+ xdx' k+1 J (fc) (x fc ; x' fc ) (5p(xj - x k+x ) - Spix'j - x k+1 )) 6 v (x k+ x - x k+1 ) 

x (^N^HxkiXk+ilx'frXk+l) -7S 1) (xfe,SA ! +i;x' A .,x' A . +1 )) -> 

for N — > oo, uniformly in s. Using (|4.55|) . (|4.57|) . and (|4.59|) . it follows from (|4.54l) . that 

k . 

(J {k \^) = (J [k) n { Z)-*lb [ ds [ dx fc dx' fc j( fe )(x fc ;x' fc )(-A i + A;.)7g s (x fc ;x' fc ) 

j=l J o J 

-ib^2 ds dx fc dx' fc dx fc+ idx' fc+1 J (fe) (x fc ; x'J^Xj - x fc+ i) - <^(x' - x fc+ i)) 
i= i-/o 7 

(4.60) 

for any fixed i and A;. Finally, we apply Lemma T4.2I to replace 5 v (x k +\ — x' k+l ) by 5(x k+ i — x' k+l ) 
and ^(xj - x fc+ i) (respectively, 5p(x' j - x k +i)) by 8(xj - x k+ i) (respectively, by 5(x' j - x k+1 )). The 
error here is of order /3 1 / 2 + rj. Hence, letting rj — > and /?^0we find 



(•^tS) = (J {k) ,l ( Z) ~*pj* ds J d ^ d *'k ^(^4) (-A, + A;.) 7 g s (x fc ;x' fc ) 
-ib^ ds dxfcdx'fcdxfc+i (x k ;x.' k )(5(xj - x k+l ) - 5(x' j - x k+1 )) 



(4.61) 



> X k +l] Xi, . . . , Xj,, Xfc+l) . 

A Sobolev Type Inequalities 

Lemma A.l. i) Suppose V G L 3 / 2 (A), and tp £ W 1,2 (A). Then 



□ 



dx|V(x)| 2 lu(x/a) < C\\V\\ L3/HA) (HWII 2 + ||V|| 2 ) 1/2 • (A.62) 
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ii) Suppose V G -^ 1 (A). Then, considering V(x — y) as an operator on L 2 (A,dy) ® L 2 (A,dx) we 
have the operator inequality 
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569-615. 



19 



